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AN  EXAMINATION  OF  THE  EFFECTS  OF  CLOUD  SEEDING 
IN  PHASE  II  OF  THE  SANTA  BARBARA  CONVECTIVE  BAND 
SEEDING  TEST  PROGRAM 

Ralph  A.  Bradley,  Thomas  C.  Redman  and  Thomas  A.  Gleeson 

Florida  State  University 
Tallahassee,  Florida 

SUMMARY 

This  report  covers  statistical  analyses  of  the  experimental  data  from  Phase  II  of 
the  Santa  Barbara  Convective  Band  Seeding  Test  Program  conducted  from  1970  to  1974. 
Comparisons  are  made  with  earlier  analyses  of  the  Phase  I  data. 

The  Phase  II  study  was  in  two  parts,  essentially  separate  experiments,  one  using 
ground-seeding  techniques  and  one  using  aerial -seeding  techniques.  Data  summaries  of 
both  precipitation  responses  and  potential  concomitant  variables  are  given  in  an 
appendix.  The  main  analyses  for  examination  of  the  effects  of  seeding  are  weighted 
analyses  of  variance  of  transformed  precipitation  data  for  various  defined  target 
areas  in  Section  5.  The  experiments  are  relatively  small  and  no  effects  of  seeding 
are  apparent  except  for  the  aerial-seeded  part  of  the  experiment,  when  border-line 
one-sided  significances  are  obtained  after  omission  of  four  storms,  not  treated 
fully  in  accordance  with  the  design  plan.  The  use  of  concomitant  variables  as 
covariates  in  covariance  analyses  is  examined,  with  tables  given  in  the  appendix. 

The  use  of  covariates  enhanced  apparent  treatment  effects  for  the  ground. seeded  part 
of  the  experiment  but  not  for  the  aerial -seeded  part. 

The  Phase  II  study  used  storms  as  the  basic  experimental  unit  whereas  the  Phase 
I  study  used  the  convective  band.  Difficulties  arise  in  analyses  because  of  this 
change.  The  covariates  were  again  measured  in  the  area  of  expected  response  from 
seeding  and  hence  are  suspect.  Improved  design  of  future  similar  studies  would  require 
use  of  better  covariates,  avoidance  of  possible  seeding  effects  on  the  covariates, 
choice  of  good  control  areas,  better  selection  of  acceptable  experimental  units, 
the  use  of  convective  bands  as  experimental  units,  if  possible  persistence  effects 
can  be  discounted,  and  the  use  of  experimental  designs  that  allow  for-  storm  effects 
as  a  source  of  variation. 
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1.  INTRODUCTION 


The  Santa  Barbara  Convective  Band  Seeding  Test  Program  was  conducted  by  North 
American  Weather  Consultants  and  its  affiliated  organization,  Aerometric  Research 
Inc.,  from  1967  through  1974.  The  test  program  was  in  two  phases,  Phase  I  involving 
the  1967-68  through  1970-71  seasons  of  winter  s'  rms,  and  Phase  II  involving  the 
1971-72  through  1973-74  seasons  with  some  preliminary  data  from  the  1970-71  season. 
Research  was  sponsored  on  two  concurrent  projects  by  the  Naval  Weapons  Center, 

China  Lake,  California  and  the  Bureau  of  Reclamation,  Denver,  Colorado,  the  second 
study  based  on  an  augmented  network  of  raingages.  Basically,  the  studies  involved 
cloud  seeding  of  convective  bands  within  winter  storms  in  the  Santa  Barbara  area 
of  California,  with  precipitation  responses  attributable  to  these  convective  bands 
recorded  by  extensive  networks  of  raingages.  Two  final  reports  were  issued,  Thompson, 
Brown  and  Elliott  (1975)  and  Brown,  Thompson  and  Elliott  (1975),  both  of  which  have 
extensive  bibliographies,  including  interim  reports  and  publications  associated  with 
the  projects. 

Research  sponsored  by  the  Office  of  Naval  Research  at  the  Florida  State  Univer¬ 
sity  was  motivated  by  the  availability  of  the  extensive  data  available  from  the  test 
program,  a  perceived  need  for  further  statistical  analysis,  and  perhaps  the  potential 
economic  importance  of  effective  cloud  seeding  of  West  Coast  winter  storms.  Methods 
of  analysis  initially  used  are  described  in  the  two  cited  final  reports  and  summar¬ 
ized  also  by  Bradley,  Srivastava  and  La.nzdorf  (1977a),  and,  in  less  detail,  by 
Bradley  and  Scott  (1T30).  The  original  statistical  analyses  were  based  on  applica¬ 
tion  of  the  Wilcoxon-Mann-Whitney,  two-sanple,  rank  test,  separately  for  the  data 
from  each  raingage.  Both  precipitation  measurements  and  ratios  of  precipitations 
at  the  raingage  divided  by  average  corresponding  control  area  precipitations  were 


used  with  similar  results  for  ground-seeded  convective  bands.  For  aerial -seeded 
bands  in  the  Phase  II  study,  no  control  area  was  used  and  analyses  were  based  on 
the  raingage  precipitations.  Extensive  new  analyses  of  the  Phase  I  data,  many  of 
then  exploratory  in  nature,  have  been  conducted  under  the  present  contract. 

Our  studies  of  the  Phase  I  data  have  been  extensive.  It  was  desired  to  develop 
a  procedure  for  data  summarization  of  precipitation  attributable  to  a  convective 
band,  the  experimental  unit  of  the  Phase  I  experiment,  for  a  designated  target  area. 

A  response  surface  method  was  used;  it  was  found  necessary  to  represent  response 
with  a  cubic  surface.  An  integrated  precipitation  response  was  developed,  but  found 
to  be  very  highly  correlated  with  the  simple  precipitation  mean,  the  latter  then 
being  considered  as  the  adequate  measure  of  precipitation  for  a  target  or  control 
area.  Efforts  on  data  summarization  have  been  reported  by  Bradley,  Srivastava  and 
Lanzdorf  (1977a,  1977b).  Gleeson  (1977)  summarized  data  on  possible  covariates 
associated  with  each  convective  band,  covariates  that  were  based  both  on  cloud  physics 
measurements  taken  by  radiosonde  and  band  passage  time  or  duration.  Bradley,  Srivastava 
and  Lanzdorf  (1979a, b)  presented  analyses  of  variance  and  covariance  of  band  precip¬ 
itation  means  and  integrated  precipitations  for  various  target  areas  using  those  co¬ 
variates  along  with  precipitation  measures  of  a  defined  control  area.  It  was  found 
that  the  standard  deviation  of  precipitations  among  raingages  in  a  target  area  was 
related  linearly  to  the  precipitation  mean;  this  suggested  analysis  of  transformed 
precipitations,  z  *  log(A*y),  where  y  is  precipitation  at  a  raingage,  z  is  the  trans¬ 
formed  response,  and  A  is  a  constant  associated  with  a  target  area,  estimated  from 
the  observed  linear  relationship.  The  transformation  was  shown  to  effectively  stabil¬ 
ize  variances,  except  perhaps  at  very  low  precipitation  levels.  Analyses  of  variance 
and  covariance  were  done  with  the  transformed  precipitations  also.  Bradley  and  Scott 
(1979,  1980),  concerned  with  the  validity  of  parametric  assumptions,  sampled  the 
randomization  distributions  associated  with  certain  of  the  analyses  and  verified  that 
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parametric  analyses  approximated  the  randomization  tests  well.  One  finding  was 
that  use  of  the  covariates  was  suspect,  the  covariates  apparently  having  been 
affected  by  treatment  (seeding);  indeed,  the  cloud  physics  covariates  were  measured 
at  Santa  Barbara  Airport,  well  into  the  intended  area  of  expected  response  to 
seeding,  and  this  oust  be  judged  a  defect  in  the  experimental  design.  Simple 
analyses  of  variance  for  defined  target  areas  yielded  one-sided  significance  levels 
of  approximately  0.06  for  the  transformed  data,  consistent  with  a  randomisation 
analysis  check  by  Elliott  and  Brown  (1971).  Scott  (1979)  used  principal  components 
for  data  summarization  with  limited  success  and  difficulties  in  application.  His 
first  component  was  highly  correlated  with  the  precipitation  mean  and,  although  two 
additional  components  could  be  identified  and  interpreted,  they  contributed  little 
in  explanation  of  variability  among  raingages  in  a  target  or  control  area. 

In  this  report,  the  more  promising  of  the  methods  of  analysis  used  for  the 
Phase  I  data  are  applied  in  analysis  of  the  Phase  II  data.  If  one  viewed  the 
Phase  II  experiment  as  a  confirmatory  one  to  verify  the  preliminary  suggestions 
of  an  effect  of  cloud  seeding  exhibited  by  the  Phase  I  experiment,  we  would  regard 
the  appropriate  analyses  to  be  those  reported  in  Table  IV  below.  But  we  report 
also  analyses  of  covariance  that  provide  conflicting  impressions  of  the  appropriate 
ness  of  the  use  of  the  selected  covariates.  There  were  design  changes  for  the  Phase 
II  experiment  and  these  design  changes  detract  from  the  experiment  as  a  confirmatory 
one.  There  are  anomalies  in  the  data  that  may  suggest  the  need  for  more  care  in 
the  acceptance  of  a  convective  band  as  an  experimental  unit. 

In  Section  2  of  this  report,  design  changes  in  the  experiment  are  highlighted 
and  discussed.  Section  3  describes  the  data  summarization  used,  while  Section  4 
explains  the  methods  of  analysis.  Analyses  of  variance  for  the  effect  of  seeding 
are  given  and  discussed  in  Section  5,  together  with  discussion  of  analyses  of  covar¬ 
iance  exhibited  in  Appendix  tables.  The  report  concludes  with  discussion  and  remarks 
on  the  experimentation  and  some  comments  on  the  improved  design  of  future  similar 
experiments. 
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2.  PHASE  II  DESIGN  CHANGES 

Phase  I  and  Phase  II  of  the  Santa  Barbara  Convective  Band  Seeding  Test  Program 
were  designed  very  similarly.  Both  phases  involved  cloud  seeding  of  convective 
bands,  measurement  of  precipitation  by  essentially  the  same  network  of  raingages, 
and  measurement  of  cloud  physics  covariates  by  radiosonde,  if  possible  at  Santa 
Barbara  Airport  at  Goleta,  California.  Raingage  locations  for  the  Naval  Weapons 
Center  study  are  given  in  Figure  2-4  by  Thompson,  Brown  and  Elliott  (1975)  and  for 
the  Bureau  of  Reclamation  Study  in  Figure  2-4  by  Brown,  Thompson  and  Elliott  (1975). 
Comparison  of  these  figures  shows  the  extended  area  and  augmented  network  of  rain- 
gages  for  the  second  study.  But  there  were  design  changes  for  the  Phase  II  study, 
some  of  them  crucial  to  appropriate  statistical  analysis.  We  note  the  most  important 
of  the  design  changes,  while  assuming  that  the  reader  is  familiar  with  the  general 
Phase  I  experiment  from  analyses  that  we  have  reported  earlier  and  may  turn  to  the 
cited  final  reports  if  detail  is  required. 

A  decision  to  investigate  the  effects  of  aerial  seeding  in  the  Phase  II  experi¬ 
ment  was  made.  The  result  was  that  the  data  for  the  Phase  II  experiment  should  be 
considered  in  two  parts,  aerial- seeded  and  ground-seeded,  and  this  is  done  in  this 
report.  Ground-seeded  data,  seeding  done  from  the  same  seeding  site  as  in  the 
Phase  I  experiment,  resulted  when  aerial  seeding  was  not  possible.  A  randomized 
decision  to  seed  or  not  seed  was  applied  and  the  ground-seeded  data  consist  of 
responses  to  both  seeded  and  not-seeded  convective  bands  as  they  did  in  the  Phase  I 
data.  The  aerial-seeded  data  are  similar  to  the  ground- seeded  data,  with  responses 
to  both  seeded  and  not-seeded  convective  bands.  The  seeding  aircraft,  flew  its  desig¬ 
nated  flight  paths  in  both  situations  but  performed  cloud-seeding  only  on  the  appro¬ 
priate  randomized  decision.  The  seeding  aircraft  flew  at  or  near  the  freezing  level 
along  a  30  to  60  km  track  within  the  convective  band  and  transverse  to  its  direction 
of  movement.  This  was  done  in  an  area  10  to  30  km  west  of  the  coast,  upwind  of  the 
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centroid  of  the  instrumented  area.  Ground-seeded  data  resulted  when  range  scheduling 
conflicts  or  other  problems  made  aerial-seeding  impractical;  the  path  of  the  seeding 
aircraft  was  in  restricted  air  space  associated  with  Vandenberg  Air  Force  Base. 

Tables  in  this  report  are  labelled  (a)  or  (b)  for  ground-seeded  and  aerial-seeded 
data  summaries  respectively. 

The  silver  iodide  seeding  generator  was  changed  for  the  Phase  II  experimenta¬ 
tion,  with  a  new  generator  developed  by  North  American  Weather  Consultants.  Some 
changes  in  concentration  of  the  Agl-NH^I-accetone  solution  were  made  for  the  1972-73 
and  1973-74  seasons.  A  ground-based  version  of  the  airborne  acetone  burner  was 
employed  when  ground  seeding  was  necessary.  During  the  1972-73  season,  two  of  the 
seeding  flights  utilized  droppable  WMU-1/8  pyrotechnic  flares.  Aerial-seeded  data 
sets  used  in  this  report  and  in  analyses  in  the  two  final  reports  of  the  experimenters 
contain  all  aerial-seeded  convective  bands.  In  the  Phase  I  experimentation  and  in 
the  ground- seeded.  Phase  II  experimentation,  precipitation  measurements  west  of  the 
ground -seeding  site  may  be  and  have  been  used  to  provide  a  covariate,  control-area 
precipitation,  in  certain  covariance  analyses.  This  is  not  possible  for  the  aerial- 
seeded  experimentation.  Indeed,  for  analyses  of  the  aerial-seeded  data,  target 
areas  may  be  defined  that  includes  raingages  regarded  to  be  in  the  control  area  for 
the  ground- seeded  data. 

The  major  design  change  affecting  appropriate  statistical  analysis  was  that 
a  storm  became  the  experimental  unit  rather  than  the  convective  band.  We  find  the 
stated  reason  somewhat  obscure  and  quote  from  Thompson,  Brown  and  Elliott  (1975): 

"One  criticism  of  the  band-by-band  randomization  scheme  employed  in  the 
preceding  phase  of  the  program  was  that,  although  it  did  provide  ade¬ 
quately  for  the  possibility  of  interactive  effects  between  a  seeded 
band  and  unseeded  bands  preceding  or  following  it  within  a  given  storm, 
it  did  not  permit  the  testing  for  any  multiplication  effects  which 
might  occur  if  all  bands  within  a  given  frontal  zone  were  seeded. 
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To  meet  this  objection,  a  randomization  mode,  based  upon  a  rigid  48-hour 
time  block,  was  adopted  in  which,  during  the  48-hour  period  subsequent  to 
the  onset  of  precipitation,  each  convective  band  was  treated  in  accordance 
with  the  randomized  decision  for  the  block  as  a  whole.  Since  storms  in 
this  area  have  typical  durations  of  between  twelve  and  thirty-six  hours, 
this  provided  effective  randomization  on  a  storm-by- storm  basis,  while 
retaining  the  advantage  of  large  sample  size  provided  by  statistical  treat¬ 
ment  of  rainfall  data  for  individual  bands." 

We  are  not  sure  of  the  interpretation  of  "multiplication  effect"  in  the  quotation. 

It  would  appear  to  relate  to  a  build-up  effect  of  seeding  as  the  storm  progresses 
and  the  several  convective  bands  pass  over  a  target  area.  But  this  implies  a 
carry-over  effect  from  seeding  one  band  to  successive  bands  in  a  storm  sequence 
of  bands.  If  this  occurs,  it  seems  likely  to  occur  as  a  carry-over  effect  on  a 
subsequent  unseeded  band  also  and  we  cannot  reconcile  this  with  the  earlier  state¬ 
ment  that  provision  was  made  for  "the  possibility  of  interactive  effects".  Bradley, 
Srivastava  and  Lanzdorf  (1979a)  looked  for  carry-over  effects  on  unseeded  bands 
following  seeded  bands  in  the  same  storm  and  also  for  effects  due  to  positions  of 
the  band  in  the  storm  without  success  -  see  Tables  SA-V  to  SA-IX.  No  analyses  are 
given  in  the  two  final  reports  of  the  experimenters  for  the  detection  of  a  multipli¬ 
cation  effect.  Indeed,  analyses  given  in  these  two  final  reports  for  the  Phase  II 
data  parallel  those  used  for  the  Phase  I  data,  are  again  on  a  raingage-by-raingage 
basis,  and  take  no  cognizance  of  the  fact  that  a  storm  has  now  become  the  experimen¬ 
tal  unit,  contrary  to  implications  of  the  last  sentence  of  the  quotation. 
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We  shall  designate  the  variability  among  convective  bands  treated  alike  within 
storms  as  a  sampling  error.  This  is  appropriate  if  there  are  no  positional  or 
multiplication  effects.  We  shall  designate  the  variability  among  storms  treated 
alike,  seeded  or  not  seeded,  as  an  experimental  error.  If  there  is  a  component  of 
variability  in  this  experimental  error  attributable  to  variation  among  storms,  we 
would  expect  that  the  experimental  error  would  be  larger  than  the  sampling  error. 

But  variability  between  treatments,  seeded  and  not  seeded,  must  contain  also  between- 
storm  variability  and  sampling  variability,  along  with  a  possible  effect  of  treat¬ 
ment.  In  an  analysis  of  variance  or  covariance  table,  mean  squares  for  sampling 
error,  experimental  error,  and  treatment  must  be  properly  calculated.  The  appro¬ 
priate  test  for  treatment  or  seeding  effect  must  take  into  account  both  wlthin-storm 
band  to  band  variability  and  between-storm  variability. 

As  in  the  Phase  I  experiment,  covariates  were  observed  for  each  convective 
band.  When  the  experimental  unit  is  the  storm,  this  produces  a  somewhat  unusual 
analysis.  The  covariates  must  be  used  to  "standardize"  the  sampling  unit  (convective 
band)  responses;  thus  it  would  be  expected  that  the  use  of  covariates  would  reduce 
sampling  error.  They  may  be  expected  also  to  affect  mean  squares  for  experimental 
error  and  treatment,  since  these  mean  squares  have  expectations  with  a  component  of 
variability  for  sampling  error.  If  covariates  are  correlated  with  treatment,  the 
problem  encountered  in  the  Phase  I  covariance  analyses,  their  use  may  again  remove 
whatever  treatment  effect  is  present. 

Our  methods  of  analysis  are  described  in  Section  4.  It  will  be  seen  that  they 
provide  means  of  calculation  of  mean  squares  for  sampling  error  and  experimental 
error.  The  raingage-by-raingage,  Wilcoxon-1Mann-’".utney  rank  tests  of  tne  experimen¬ 
ters '  reports  are  designed  to  detect  location  changes  in  the  distributions  of  seeded 
and  non-seeded  responses  and  assume,  under  the  hypothesis  tested,  that  seeded  and 
non-seeded  responses  come  from  a  single  response  distribution.  This  is  not  the 
case  if  storm  effects  are  present. 
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3.  PHASE  II  DATA  SUMMARIZATION 

Throughout  Phase  II  data  summarization  and  analysis,  the  ground- seeded  and 
aerial-seeded  data  are  considered  separately.  Given  that  storms  are- experhuuirtal 
units,  the  tabulation  of  storms  seeded  and  not  seeded  is  as  follows: 


No.  of  Storms 

Seeded 

Not  Seeded 

Total 

Gd.  Seeded  Part 

7 

S 

12 

Aerial  Seeded  Part 

12 

18 

30 

Similarly,  on  a  convective  band  basis,  the  tabulation  is: 


No.  of  Bands 

Seeded 

Not  Seeded 

Total 

Gd.  Seeded  Part 

20 

10 

30 

Aerial  Seeded  Part 

18 

27 

45 

Some  notes  are  in  order.  A  storm  and  its  convective  bands  became  part  of  the 
ground-seeded  or  aerial-seeded  experimentation  because  of  the  impractical ity  or 
practicality  of  aerial  seeding;  the  randomized  decision  on  whether  or  not  to  seed 
came  later.  Thus  the  ground-seeded  part  of  the  experimentation  included  both  seeded 
and  not-seeded  storms  as  did  the  aerial- seeded  part  of  the  experiment.  In  addition, 
there  were  four  storms  that  were  part  of  the  aerial -seeded  data  with  some  bands 
seeded  and  some  bands  not-seeded  -  see  Table  A. 1(b),  storms  91,  92,  3  and  12.  In 
the  table  above  and  in  some  analyses,  each  of  these  storms  was  treated  as  two  storms, 
one-seeded  and  one  not-seeded. 

The  precipitation  data  are  summarized  in  Tables  A- 1(a)  and  A-l(b)  of  the  appendix 


by  storms,  bands  and  response  areas.  For  each  convective  band  in  each  storm  for 
each  response  area,  the  mean  precipitation  in  inches  is  given  along  with  the  number 
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of  operative  raingages  in  the  response  area  and  the  variance  among  those  raingages. 
The  response  areas  are  defined  in  Table  I,  the  first  five  of  which  were  used  also 
in  our  Phase  I  data  analyses.  Response  area  (i)  is  the  main  target  area  for  ground- 
seeded  experimentation.  Areas  (ii)  and  (iii)  are  respectively  near  and  far  from  the 
ground  seeding  site  and  are  used  to  investigate  areas  of  effect  of  ground  seeding. 
Area  (v)  is  a  control  area  west  and  up-wind  of  the  seeding  site  for  ground-seeded 
experimentation;  but  becomes  a  near  to  seeding  target  area  for  aerial  seeding. 

Area  (vi)  is  a  total  target  area  for  aerial  seeding.  Data  summaries  for  response 
areas  (i)  -  (vi)  are  based  on  the  augmented  network  of  raingages  used  in  the  Bureau 
of  Reclamation  study.  Areas  (vii)  and  (viii)  are  control  and  target  areas  for  the 
Naval  Weapons  Center  Study  for  ground-seeded  data,  while  both,  along  with  Area 
(ix)  are  target  areas  for  aerial-seeded  data.  Data  summaries  and  analyses  for 
Areas  (vii)  -  (ix)  are  based  on  the  raingage  network  of  the  Naval  Weapons  Center 
Study.  The  numbers  of  stations  in  Table  I  are  the  numbers  of  raingages  in  the 
response  areas,  not  all  of  which  were  always  operative.  Reference  to  Figures  2-4 
of  the  two  final  reports  will  assist  in  understanding  the  defined  response, 
areas.  Tables  A-l(a)  and  A-l(b)  do  not  have  data  for  response  area  (ix),  but  such 
data  may  be  reconstructed  from  those  for  Areas  (vii)  and  (viii)  if  desired. 
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TABLE  I 

Definitions  of  Response  Areas* 


Response 

Area 

Ranges 

Latitude 

in  Degrees 

Longitude 

Number  of 
Stations 

(i) 

34.0  -  35.25 

118.0  -  120.02 

106 

(ii) 

34.4  -  35.0 

119.51  -  120.02 

25 

(iii) 

34.0  -  35.0 

118.0  -  119.51 

71 

(iv) 

Areas  (ii)  +  (iii) 

96 

(v) 

34.4  -  35.25 

120.02  -  120.60 

34 

(vi) 

Areas  (i)  ♦  (v) 

140 

(vii) 

All  stations  in  Naval  Weapons  Center  Study 

41 

West  of  seeding 

site  at  120.02°  long. 

(viii) 

All  stations  in 

Naval  Weapons  Center  Study 

63 

East  of  seeding  site  at  120.02°  long. 

(ix) 

Stations  and  Areas  of  (vii)  +  (viii) 

104 

Precipitation  means  in  inches  for  the  various  response  areas  are  given  in 
Tables  11(a)  and  11(b)  respectively  for  ground-seeded  and  aerial-seeded  bands. 

We  have  exhibited  these  means  for  visual  comparisons  of  seeded  and  not-seeded 
results,  since  our  analyses  are  based  on  transformed  precipitations  for  which  means 
are  less  easily  interpreted. 

Examination  of  Table  A-l(a)  shows  that  not-seeded  storm  4  with  one  band  has 
extremely  high  precipitation.  We  have  treated  it  as  an  outlier  and  performed 
analyses  with  this  storm  included  and  excluded.  The  first  section  of  Table  11(a) 
gives  means  for  all  ground-seeded  bands  while  the  second  section  omits  this  storm 
and  band.  It  is  seen  that  the  omission  of  the  extreme  storm  has  a  major  effect  on 
the  apparent  effect  of  seeding.  Indeed,  further  examination  of  Table  A- 1(a)  sug¬ 
gests  that  non-seeded  band  9-1  has  high  precipitation;  while  seeded  bands  8-2  and 
22-3  have  somewhat  high  precipitations.  These  notes  reflect  the  high  variability 
of  cloud  seeding  data;  in  this  small  ground-seeded  experiment, the  'luck  of  the  draw" 
in  the  randomization  may  play  an  important  role.  From  another  viewpoint,  this  vari¬ 
ation  may  suggest  need  for  better  criteria  for  determination  of  a  "seedable"  band. 
Table  11(b)  is  in  two  parts,  the  first  with  all  storms  and  bands  and  the  second 
with  the  four  mixed  storms  discussed  above  omitted.  It  happened  that,  for  all 
four  storms,  the  unseeded  bands  had  higher  precipitation  means  than  the  seeded  bands. 
We  do  not  like  the  removal  of  data  without  assignable  cause  and  tend  to  put  most 
weight  on  analyses  with  all  of  the  available  data  used.  The  means  in  Tables  11(a) 
and  11(b)  are  simple  means  of  the  convective  band  means  of  Tables  A-l(a)  or  Tables 
A- 1(b)  as  appropriate. 


I 
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Four  storms  had  both  seeded  and  unseeded  bands.  In  Part  1,  seeded  and  unseeded  bands  in  sequences  were  treated 
as  different  storms;  in  Part  2,  these  storms  were  omitted.  See  text. 


For  analyses  of  variance  and  covariance,  as  in  the  Phase  I  analyses  of  Bradley, 
Srivastava  and  Lanzdorf  (1979a,  b),  it  seems  appropriate  to  transform  the  data  to 
stabilize  variances.  This  was  done  by  plotting  the  standard  deviation  among  rain- 
gage  observations  for  a  band  and  response  area  against  the  corresponding  mean  res¬ 
ponse,  the  data  being  available  from  Tables  A-l(a)  and  A-l(b).  It  is  found  that 
these  plots  are  very  nearly  linear.  Straight  lines  were  fitted  by  least  squares 
and  slopes  and  intercepts,  B  and  A,  were  estimated.  The  transformation  was 
z  =»  log(A+y)  where  y  is  precipitation  at  a  raingage  and  A  *  A/B.  This  procedure 
was  illustrated  for  Phase  I  data  in  Figures  4  and  5  of  Bradley,  Srivastava  and 
Lanzdorf  (1979b) .  Figures  for  the  various  target  areas  for  ground-seeded  and 
aerial-seeded  data  for  Phase  II  were  very  similar.  The  transformation  does  stabil¬ 
ize  variances,  except  perhaps  for  very  small  values  of  y.  Values  of  A  are  given 
in  Tables  111(a)  and  111(b).  Ground-seeded  values  in  Table  111(a)  are  slightly 
larger  than  those  found  in  the  Phase  I  analyses,  the  latter  ranging  from  0.03  to 
0.075.  The  basic  data  summarization  of  precipitation  data  in  preparation  for 
analyses  of  variance  and  covariance  involves  preparation  of  tables  like  Tables  A-l(a) 
and  A-l(b)  giving  means  of  the  transformed  raingage  precipitations  by  storms  and 
bands  for  the  various  response  areas.  It  was  a  secondary  effect  of  the  use  of  the 
transformation  that  the  influence  of  larger  precipitation  bands  was  somewhat 
reduced;  weighted  means  of  transformed  data  were  such  that  the  means  of  seeded  bands 
were  larger  than  means  of  not-seeded  bands  for  all  but  one  of  the  target  areas  - 
see  the  signs  in  Tables  IV(a)  and  IV(b)  in  Section  5. 
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Gleeson  (1977)  defined  and  summarized  data  on  concomitant  variables  or  covari¬ 
ates  for  the  Phase  I  data.  Similar  data  summaries  are  provided  in  Tables  A-2(a) 
and  A-2(b)  of  the  appendix  of  this  report  for  the  Phase  II  data.  Table  notes  on 
these  tables  are  given  also  in  the  Appendix,  providing  explanations  of  these  concom¬ 
itant  variables.  Covariates  considered  are  listed  in  order  below  with  brief  descrip 
tive  labels,  their  designations  matching  those  used  in  the  Phase  I  reports: 

X^t  Mixing  Ratio, 

X2:  700  mb  Wind  Speed, 

X3:  700  mb  Wind  Direction, 

X^ :  Mean  Wind  Speed, 

X,.:  Direction,  Avg.  Vector  Wind, 

X^:  500  mb  Temperature,  (3.1) 

X7:  Stability  Class, 

Xg:  Showalter  Index, 

Xg:  Stability  Wind  Speed, 

X1Q:  Direction,  Stability  Wind, 

X^:  Instability  Transport, 

X^:  Band  Passage  Time  (Seeding  Site),  Duration. 

The  data  in  Tables  A-2(a)  and  A-2(b)  are  used  with  the  band  means  of  the  transformed 
precipitation  data  for  response  area  (i)  to  exhibit  correlations  in  Tables  A-3(a) 
and  A-3(b)  for  Area  (i).  In  Table  A-3(a),  X£  and  Z  appear  also.  X£  is  the  pre¬ 
cipitation  mean  by  bands  for  response  area  (v) ,  the  control  area  for  ground  seeding, 
the  means  based  on  the  Bureau  of  Reclamation  data.  Z  is  the  seeding  indicator  vari¬ 
able.  These  tables  are  given  for  comparison  with  Table  A-l  of  Bradley,  Srivastava 


and  Lanzdorf  (1979a).  The  patterns  of  correlations  are  very  similar,  perhaps 
somewhat  larger  in  magnitude  for  Phase  II  data.  Two  correlations  merit  comment. 

The  correlations  between  transformed  precipitation  band  mean  and  Xc>  control  area 
precipitation,  are  larger,  changing  from  0.7  to  0.9.  The  correlations  between 
transformed  precipitation  banc!  mean  and  X^2>  band  passage  time  or  duration,  are 
smaller,  particularly  for  the  aerial -seeded  data. 

Detailed  analyses  of  the  Phase  I  data  showed  redundancies  in  the  covariates. 

It  was  judged  that  the  set  of  covariates  could  be  reduced  to  X2,  X^,  Xft,  X?,  Xg, 
and  Xjj  along  with  consideration  of  Xc  and  X^2>  This  is  done  again  in  this  report 
in  the  covariance  analyses.  We  shall  regard  the  first  six  covariates  as  the  basic 
set.  Models  considered  in  analyses  of  variance  and  covariance,  differently  labelled 
than  in  Phase  I  analyses,  are  as  follows: 

Model _ Identification  of  Covariates _ 

(1)  No  covariates 

Uj  X2’  X3’  X6’  X7*  V  X11 

(3)  Covariates  of  (2)  plus  X^2  (3.2) 

(4)  Covariates  of  (2)  plus  X£. 

(5)  Covariates  of  (2)  plus  X^2,  Xc> 

Description  of  data  summarization  for  analyses  reported  in  subsequent  sections  of 
this  report  is  now  complete.  The  methods  of  analysis  are  described  in  the  follow¬ 
ing  section. 
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4.  METHODS  OF  STATISTICAL  ANALYSIS 

Bradley,  Srivastava  and  LanzdoTf  (1979a, b)  developed  analyses  of  variance 
and  covariance  by  target  areas  for  transformed  precipitations  and  covariate 
values  using  the  convective  band  as  the  experimental  unit.  Since  the  storm  has 
become  the  experimental  unit  in  the  Phase  II  experimentation,  changes  in  the 
methods  of  analysis  are  necessary  and  our  new  methods  are  explained  below. 

Consider  first  the  analysis  of  variance.  The  basic  linear  model  is 


z.  .  =  p  ♦  t.  +  3. e.  .  , 

13a  i  3(1)  13a' 


(4.1) 


where  z^ja  is  the  mean  of  the  transformed  precipitations  at  raingages  in  a  target 

area  for  band  a  of  storm  j  receiving  treatment  i,  o  a  l,...,n.^,  n^  being  the 

number  of  bands  in  storm  j  with  treatment  i,  3  »  l,...,m,  n^  being  the  number  of 

storms  receiving  treatment  i,  and  i  =  1,2,  for  seeded  and  not  seeded  storms  respec 

tively.  The  components  of  the  model  in  the  right-hand  side  of  (4.1)  are  y,  the 

general  mean,  the  effect  of  treatment  i,  the  effect  of  storm  3'  of  treat 

ment  i  taken  as  a  member  of  a  set  of  independent  random  variables  with  zero  means 
2 

and  variances  a  ,  and  a  member  of  a  second  set  of  independent  random  vari- 

2 

ables,  independent  of  the  first  set,  with  zero  means  and  variances  where 

nn  is  the  number  of  raingages  contributing  to  z„o. 

Transformation  of  raingage  precipitations  was  effected  to  stabilize  variances 

2 

permitting  the  assumption  that  the  variance  of  e. .  is  proportional  to  a  / m, .  , 

X  J  01  S  X  J  01 

since  1jL^a  *s  the  mean  of  nu ^  transformed  raingage  precipitations.  The  variance 
2 

of  would  be  CTs/mjja»  if  each  transformed  raingage  precipitation  had  variance 


■* 

1 


18 


og  and  such  precipitations  in  a  target  area  for  a  convective  band  were  independent; 
we  believe  that  the  proportionality  assumption's  adequate.  To  develop  the  analysis 
of  variance  table  for  the  model  (4.1),  it  is  necessary  to  use  weighted  least  squares 
in  minimization  of 


l  l  fjm..  (z..  -  p  -  T.  - 

>1  ^  13a  1  j(i); 


(4.2) 


There  is  a  redundancy  of  treatment  parameters  and  it  is  convenient  to  add  the  con- 
2 

straint  that  £  m.  r  =  0  and  to  minimize  (4.2)  subject  to  this  constraint.  The 
i-1  1 

number  of  raingage  observations  contributing  to  estimation  of  the  mean  of  the  trans- 
ftarmSd  preciprtatiofe  for  treatment  i  is  ra  .  The  analysis  of  variance  table  result¬ 
ing  is  as  follows: 


Source  of 
Variation 

d.f. 

s.s. 

Expected  M.S. 

Treatment 

1 

2  .  7 

I  m  (z  -i  y 

-1  *  •  A  *  *  •  •  • 

2  .  2  r  2 

<T«;+k7a  +  l  mi  Ti 
i=l  * 

Experimen¬ 
tal  Error 

2 

I 

i-1 

ni 

7  m. .  (z.  .  -z  )^ 
i3.v  13.  ...' 

2  .  2 

Vkl° 

Subtotal 

w1 

2 

I 

i=l 

ni 

I  mij  (5t1  -l  )2 

j-l  13  *  13 '  •** 

(4.3) 

Sampling 

Error 

2  nj 

l  U  *»  -1) 

i-1  j=l  1J 

2 

I 

i-1 

ni  nii 

y  y  m. .  (z. .  -z. , 

j-1  a-1  130  13a  ^ 

,2  2 
i.5 

Total 

2  nj 

I  I  n..  -  1 
i-1  3-1  lj 

2 

I 

i-1 

ni  nii 

y  y  m..  (z..  -z 

j*l  a-1  • 

)2 

•  • 

Total 


Some  symbols  in  (4.3)  require  definition: 


.  ,  m. 

ii  1-- 


5i  "ij 

j=l  a=l  mija>  i3 


2  ni  n. . 

-  I  l  V\Sa  ■ 

i=l  j  =1  a=l 
nij 

=  T  Jm. .  z .  .  /m.  .  , 
0“1  iJ- 


ni  nij 

z .  =  T  V  Jm. .  z. .  /m.  ,  z 

x”  j*i  a=i  13a  *3«  I--’  ... 


2  n^  n,  • 

y  y  y  3m. .  z. .  /m  , 

i*l  j=l  o*l  tJal3a  * 


ki  s  (m...  “  £  iT2-  .^ij.^W25' 

1=1  i..  j=l  J 


(4.4) 


2  i—  — ,  n  • 

*  i  j — iq  j1.?. 

i=ib*. 


An  immediate  problem  is  seen  with  (4.3)  in  view  of  (4.4).  To  test  for  treat¬ 
ment  or  seeding  effect,  we  need  *  k^.  If  the  experiment  were  balanced  so  that 

*  m,  for  all  i,  j  and  «,  and  =  n  for  all  i,j,  then  kj  »  k2  *  mn.  If  the 

.  2 
assumption  in  the  model  were  changed  so  that  it  is  assumed  that  V'(6j  ^)  =  a  /m^  , 

then  kj  ■  k2  ■  1.  This  change  does  not  seem  appropriate. 

The  problem  is  that,  because  of  the  lack  of  balance  in  the  experiment,  the 

effects  of  storms  and  treatment  are  somewhat  confounded.  We  have  demonstrated  a 

defect  in  the  design  of  the  Phase  II  experiment,  one  that  is  implicit  in  any  analysis 

including  the  raingage-by-raingage  analyses  of  the  investigators'  final  reports. 

It  remains  to  investigate  the  extent  of  the  difficulty  and  this  is  done  in  the 

following  section. 
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Covariate  analyses  are  given  in  Tables  A-4(a)  and  A-4(b)  of  the  appendix  for 
the  ground-seeded  and  aerial-seeded  data  respectively.  General  linear  regression 
methods  were  used  to  develop  the  analysis  of  covariance  tables.  Model  (4.1)  was 
modified  to  include  covariates  X, ,...,X  kith  values  x...  ,...,x  ..  for  band  o 

*  p  iija  pija 

of  storm  j  with  treatment  i.  It  was  necessary  to  introduce  design  parameters 


also  to  represent  the  nj+n2-2  parameters  for  storm  contrasts  within  treatments  and 
these  may  be  designated  with  value  w^j  for  a  band  in  storm  j  and  treatment  i, 
k  *  l,...(n^-l).  Values  of  w^j  were  taken  proportional  the  elements  of  a  gen¬ 


eralized  Helmert  matrix: 


of  (4.1). 
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Weighted  regression  is  employed  and  the  sum  of  squares  to  be  minimized  is 


I  I  I  'in. .  (z. . 
i-1  j-1  o«l 


ni; 


x  i  p 

-u-tZj  -  X  ekw  k  -  X  r 
1  k»i  K  r~ ' 


»  X  . .  )' 
J  Y  Yija 


(4.6) 


The  analysis  of  covariance  table  is  developed  through  introducing  terns  in  the 
model  (4.5)  in  the  appropriate  sequence  and  use  of  the  corresponding  minimum  sums 
of  squares  from  (4.6).  Let  the  reduced  models  building  to  (4.5)  depend  on  subsets 
of  the  model  parameters  as  shown  in  (4.7)  with  the  terms  in  the  right-hand  side 
of  (4.5)  numbered  in  order,  and  let  the  corresponding  minimum  weighted  sums  of 
squares  be  also  designated  as  indicated  in  the  final  column  of  (4.7): 


Model 

Terms 

Sums  of  Squares 

\ 

1 

SS(A) 

B 

1.4 

SS(B) 

C 

1,2, 3, 4 

SS(C) 

(4.7) 

D 

1,3,4 

SS(D) 

The  analysis  of  covariance  table  is  obtained  from  the  indicated  sums  of  squares: 


Source  of 

Variation 

d.f. 

S.S. 

Treatment 

1 

SS(D)-SS(C) 

Experimental 

Error 

n^n^ 

SS(B)-SS(D) 

Subtotal 

nj+n2-l 

SS(B) -SS(C) 

(4.8) 

Covariates 

P 

SS (A) -SS(B) 

Sampling 

Error 

“2 

l  l  (»»-« 

i-1  j-1  13  J 

-p  SS(C) 

2  n, 

r  r* 


Total 


-1 


SS(A) 
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The  analysis  of  covariance  is  subject  to  the  same  confounding  of  storm  effects 
and  treatment  effects  as  discussed  above  in  regard  to  the  analysis  of  variance. 

If  the  fourth  term  in  the  right-hand  side  of  (4.5)  is  omitted,  then  the 
analysis  of  variance  of  (4.3)  results  from  the  weighted  regression  analysis. 

5.  ANALYSES  OF  VARIANCE  AND  COVARIANCE 

Response  areas  (i)  -  (iv)  and  (viii)  are  target  areas  for  the  ground-seeded 
experimentation  and  areas  (i)  -  (ix)  are  possible  target  areas  for  the  aerial-seeded 
experimentation.  Target  areas  (i)  and  (vi)  represent  total  target  areas  for  ground- 
and-aerial-seeded  experimentation  with  the  Bureau  of  Reclamation  data  and  target 
areas  (viii)  and  (ix)  respectively  for  the  Naval  Weapons  Center  data.  For  the 
ground-seeded  experimentation,  response  areas  (v)  and  (vii)  represent  control  areas 
and  yield  values  of  the  covariate  Xc,  control  area  mean  precipitation,  for  the 
Bureau  of  Reclamation  data  and  the  Naval  Weapons  Center  data  respectively.  There 
are  no  control  areas  for  the  aerial-seeded  experimentation. 

Various  models  for  analyses  are  given  in  (3.2).  Model  (1)  corresponds  with 
the  model  for  the  analysis  of  variance,  (4.1).  Models  (2)  -  (5)  yield  analyses  of 
covariance,  see  (4.5),  with  p  =  6,  7,  /,  and  8  respectively.  Table  IV(a)  gives 
the  analysis  of  variance  tables  for  target  areas  (i)  -  (iv)  and  (viii)  based  on 
band  means  of  transformed  target-area,  raingage  precipitations  for  all  30  bands 
and  12  storms  of  the  ground- seeded  experimentation.  Table  IV(a)  Continued  matches 
Table  IV(a)  except  that  Storm  4  with  one-band,  not-seeded  has  been  omitted  and  the 
analyses  are  based  on  29  bands  and  11  storms  -  see  Table  11(a)  and  comments  in 
Section  3.  Table  IV(b),  Part  1,  give  the  analysis  of  variance  tables  for  target 
areas  (i)  -  (ix)  based  on  band  means  of  transformed  target-area,  raingage  precipi¬ 
tations  for  all  45  bands  and  30  storms  of  the  aerial -seeded  experimentation. 
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TABLE  IV (a) 


Analyses  of  Variance  (Ho 
Che  Various  Target  Areas, 
Ground -Seeded 

Covariates) 

Transformed 

Bands 

for 

Data 

Target 

Area 

Source  of 
Variation 

d.f. 

Mean 

Squares 

F  - 

Ratio 

Sign 

(i) 

Seeding 

1 

0.14 

0.00 

♦ 

Exp.  Error 

10 

45.74 

- 

Subtotal 

11 

41.59 

- 

Sampl.  Error 

18 

15.44 

- 

(ii) 

Seeding 

1 

2.08 

0.10 

♦ 

Exp.  Error 

10 

21.00 

- 

Subtotal 

11 

19.28 

- 

Sampl.  Error 

18 

5.70 

- 

(iii) 

Seeding 

1 

0.01 

0.00 

Exp.  Error 

10 

24.74 

- 

Subtotal 

11 

22.50 

- 

Sampl.  Error 

18 

10.51 

- 

(iv) 

Seeding 

1 

0.12 

0.00 

+ 

Exp .  Error 

10 

42.81 

- 

Subtotal 

11 

38.93 

- 

Sampl.  Error 

18 

15.69 

- 

(viii) 

Seeding 

1 

1.12 

0.03 

•f 

Exp.  Error 

10 

38.42 

- 

Subtotal 

11 

35.03 

- 

Sampl.  Error 

18 

8.19 

- 

TABLE  IV(a)  -  Continued 

Analyses  of  Variance  (No  Covariates)  for 
the  Various  Target  Areas,  Transformed  Data 
Ground-Seeded  Bands,  Extreme  Stor..i  Omitted* 


Target 

Area 

Source  of 
Variation 

d.f. 

Mean 

Squares 

F  - 
Ratio 

Sign 

(i) 

Seeding 

1 

20.73 

0.71 

♦ 

Exp .  Error 

9 

29.16 

- 

Subtotal 

10 

28.32 

- 

Saaspl .  Error 

18 

16.96 

- 

(ii) 

Seeding 

1 

11.93 

0.88 

♦ 

Exp.  Error 

9 

13.59 

- 

Subtotal 

10 

13.43 

- 

Sampl .  Error 

18 

5.35 

- 

(iii) 

Seeding 

1 

9.85 

0.61 

♦ 

Exp.  Error 

9 

16.16 

- 

Subtotal 

10 

15. :3 

- 

Sampl.  Error 

18 

12.28 

- 

(iv) 

Seeding 

1 

19.47 

0.70 

♦ 

Exp.  Error 

9 

27.81 

- 

Subtotal 

10 

26.98 

- 

Sa;ap.'. .  Error 

18 

17.48 

- 

(viii) 

Seeding 

1 

21.11 

0.82 

♦ 

Exp.  Error 

9 

25.77 

- 

Subtotal 

10 

25.30 

- 

Sampl.  Error 

18 

8.58 

- 

*One  extreme  precipitation  storm,  Storm  4  with  cne  band,  not  seeded,  has  been 
omitted. 
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TABLE  IV(b) 

Analyses  of  Variance  (No  Covariates)  for 
the  Various  Target  Areas,  Transformed  Data 
Aerial-Seeded  Bands,  Part  1-All  Bands 


Target  Source  of  d.f.  Mean  F-  Sign 

Area  Variation  Squares  Ratio 


(i)  Seeding  1 

Exp.  Error  28 

Subtotal  29 

Sampl.  Error  15 

(ii)  Seeding  1 

Exp.  Error  28 

Subtotal  29 

Sampl.  Error  15 

(iii)  Seeding  1 

Exp.  Error  28 

Subtotal  29 

Sampl .  Error  15 

(iv)  Seeding  1 

Exp.  Error  28 

Subtotal  29 

Sampl.  Error  15 

(v)  Seeding  1 

Exp.  Error  28 

Subtotal  29 

Sampl.  Error  15 


80.22  0.89  ♦ 

89.74 

89.42 

80.82 

31.35  1.27  + 

24.78 

25.01 

26.44 

40.29  0.71  «• 

57.05 

56.48 

45.93 

76.16  0.90  ♦ 

84.91 

84.61 

75.11 

23.40  1.51  ♦ 

15.50 

15.77 

15.55 
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TABLE  IV(b)  -  Continued 

Analyses  of  Variance  (No  Covariates)  for 
the  Various  Target  Areas,  Transformed  Data 
Aerial -Seeded  Bands,  Part  1-All  Bands 


Target 

Area 

Source  of 
Variation 

d.f. 

Mean 

Squares 

F- 

Ratio 

Sign 

(vi) 

Seeding 

1 

87.09 

0.98 

+ 

Exp .  Error 

23 

88.85 

- 

Subtotal 

29 

88.79 

- 

Samp! .  Error 

15 

82.76 

- 

(vii) 

Seeding 

1 

18.22 

1.10 

Exp.  Error 

28 

16.49 

- 

Subtotal 

29 

16.55 

- 

Sampl.  Error 

15 

14.53 

- 

(viii) 

Seeding 

1 

80.86 

1.38 

+ 

Exp.  Error 

28 

58.77 

- 

Subtotal 

29 

59.53 

- 

Sampl .  Error 

15 

63.50 

- 

(ix) 

Seeding 

1 

62.53 

1.13 

♦ 

Exp .  Error 

28 

55.32 

- 

Subtotal 

29 

55.60 

- 

Sampl.  Error 

15 

55.19 

- 
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TABLE  IV  (b)  -  Continued 


Analyses  of  Variance  (No  Covariates)  for 
the  Various  Target  Areas,  Transformed  Data 
Aerial -Seeded  Bands,  Part  2 -Mixed  Storms  Omitted 


Target 

Area 

Source  of 
Variation 

d.f. 

Mean 

Squares 

F  - 

Ratio 

Sign 

(i) 

Seeding 

1 

320.97 

3.88* 

♦ 

Exp.  Error 

20 

82.73 

- 

Subtotal 

21 

94.07 

- 

Sampl .  Error 

13 

89.57 

- 

(ii) 

Seeding 

1 

91.18 

4.06* 

•f 

Exp .  Error 

20 

22.45 

- 

Subtotal 

21 

25.73 

- 

Sampl.  Error 

13 

29.21 

- 

(iii) 

Seeding 

1 

186.67 

3.49* 

♦ 

Exp .  Error 

20 

53.56 

- 

Subtotal 

21 

59.86 

- 

Sampl.  Error 

13 

51.10 

- 

(iv) 

Seeding 

1 

301.78 

3.86* 

+ 

Exp .  Error 

20 

78.18 

- 

Subtotal 

21 

88.82 

- 

Sampl.  Error 

13 

83.21 

- 

(v) 

Seeding 

1 

63.47 

4.44* 

♦ 

Exp.  Error 

20 

14.31 

- 

Subtotal 

21 

16.65 

- 

Sampl.  Error 

13 

17.87 

♦Values  significant  for  one- 

■sided  test 

(positive  difference 

;n  means) 

at  .OS 

level  of  significance. 


28 


TABLE  IV(b)  -  Continued 


Analyses  of  Variance  (No  Covariates)  for 
the  Various  Target  Areas,  Transformed  Data 
Aerial-Seeded  Bands,  Part  2-Mixed  Storms  Omitted 


Target 

Area 

Source  of 
Variation 

d.f. 

Mean 

Squares 

F  - 

Ratio 

Sign 

(vi) 

Seeding 

1 

327.35] 

4.02* 

+ 

Exp.  Error 

20 

81.46 

- 

Subtotal 

21 

93.17 

- 

Sampl.  Error 

13 

93.63 

- 

(vii) 

Seeding 

1 

54.98 

3.56* 

♦ 

Exp.  Error 

20 

15.43 

- 

Subtotal 

21 

17.32 

- 

Sampl.  Error 

13 

17.10 

- 

(viii) 

Seeding 

1 

232.07 

4.25* 

♦ 

Exp .  Error 

20 

54.64 

- 

Subtotal 

21 

63.09 

- 

Sampl.  Error 

13 

69.69 

- 

(ix) 

Seeding 

1 

196.49 

3.84* 

♦ 

Exp .  Error 

20 

51.19 

- 

Subtotal 

21 

58.10 

- 

Sampl.  Error 

13 

61.86 

- 

•Values  significant  for  one-sided  test  (positive  difference  in  means)  at  .05 
level  of  significance. 
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Table  IV(b),  Part  2,  matches  Table  IV (b),  Part  1,  except  that  four  storms,  sub¬ 
divided  into  separate  seeded  and  not- seeded  storms  as  described  in  Section  3,  have 
been  omitted.  The  final  column  in  Tables  IV (a)  and  IV (b)  shows  the  sign  of  the 
difference  in  the  mean  of  the  transformed  precipitations  for  seeded  bands  and  not- 
seeded  bands;  it  is  interesting  that  the  effect  of  transformation  was  to  reverse 
the  apparent  negative  effect  of  seeding  exhibited  in  Table  11(a),  except  for 
Target  Area  (iii),  when  all  30  bands  and  12  storms  are  considered.  The  transfor¬ 
mations  reduce  the  contributions  of  extreme  observations. 

Analyses  of  covariance  are  given  in  Tables  A-4(a),  A-4(a)  Continued,  A-4(b), 
Part  1,  and  A-4(b),  Part  2  for  models  (2)  -  (5),  organized  in  parallel  with  the 
analyses  of  variance  discussed  above.  No  column  is  given  for  "Sign"  since  the 
sign  is  positive,  after  adjustment  for  covariates,  in  every  analysis. 

From  experience  with  the  Phase  I  analyses,  the  various  analyses  of  variance 
were  planned  as  the  main  analyses.  The  covariance  analyses  were  done  to  check 
the  Phase  I  covariance  analyses,  when  use  of  the  covariates  was  not  helpful  and  two 
possible  covariates,  band  passage  time  or  duration  and  control  area  precipitation, 
removed  any  possible  existing  effects  of  seeding.  In  each  of  our  tables,  analyses 
of  variance  and  covariance  are  highly  correlated  since  the  same  convective  bands 
are  involved  for  all  ground-seeded  data  and  for  all  aerial-seeded  data.  The  vari¬ 
ous  target  areas  were  used  to  compare  near  and  far-away  effects  of  seeding;  we 
see  no  differences  from  one  target  area  to  another.  The  analyses  of  variance  for 
ground-seeded  data  give  no  indication  of  a  seeding  effect,  even  with  removal  of 
the  one  extreme-precipitation,  not-seeded  storm.  There  is  little  indication  of 
a  seeding  effect  when  all  aerial -seeded  bands  are  considered,  but  removal  of  the 
four  storms,  for  which  some  bands  were  seeded  and  some  not  seeded,  yielded  signifi¬ 
cances  at  the  one-sided  0.05  level  in  Table  IX(b),  Part  2,  significances  comparable 
to  those  observed  with  the  ground -seeded,  Phase  I  analyses.  While  the  four  storms 
were  removed  because  of  defect  in  following  the  experimental  plan  to  seed  or  not 


seed  all  bands  in  a  storm,  it  happened  that  all  non-seeded  bands  removed  had  higher 
mean  precipitations  than  the  seeded  bands  removed.  It  is  apparent  that  a  few  bands 
or  storms  can  have  a  major  effect  on  the  outcome  of  the  experiment.  We  believe 
that  better  control  in  the  determination  of  seedable  convective  bands  is  needed. 

It  was  noted  in  Section  4  that  there  was  some  confounding  of  treatment  and 
storm  effects  because  of  disproportionate  numbers  of  bands  in  storms  and  varia¬ 
tions  in  the  numbers  of  operative  raingages.  The  extent  of  the  problem  depends 
on  the  relative  sizes  of  k^  and  k2  in  (4.4).  We  have  calculated  k^  =  217.03  and 

k^  ■  244.33  for  Target  Area  (i),  ground-seeded  bands,  in  Table  IV(a).  In  addition, 

2  2 

considering  the  expected  mean  squares  in  (4.3),  we  estimate  os  and  a  respectively 

as  15.44  and  0.121,  The  experimental  error  mean  square  underestimates  the  error 

component  in  the  treatment  mean  square  by  3.25.  The  test  of  significance  is  biased 

in  the  direction  of  finding  an  apparent  treatment  effect.  Even  so,  significances 

were  not  found  except  in  Table  IV(b) ,  Part  2.  Let  us  similarly  examine  the  Target 

Area  (i)  analysis  of  variance  in  Table  IV(b),  Part  2.  Now  k^  *  140.05  and 

k2  s  165.228,  but  in  this  analysis  of  variance  the  experimental  error  is  smaller 

than  the  sampling  error  as  estimated  by  their  mean  squares  and  we  must  estimate 
2 

a  as  zero.  The  analyses  of  variance  in  Table  IV(b),  Part  2,  seem  appropriate 
if  deletion  of  the  mixed  storms  is  appropriate. 

A  difference  between  the  analyses  of  variance  for  ground-seeded  analyses  and 
aerial -seeded  analyses  is  now  apparent.  For  ground- seeded  data,  the  experimental 
error  is  two  to  three  times  as  large  as  the  sampling  error  and  this  was  expected 
intuitively  in  the  thought  that  variability  among  storms  treated  alike  should  be 
greater  than  variability  among  bands  treated  alike  within  the  same  storms.  But 
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this  is  not  the  case  with  the  aerial-seeded  data,  as  mean  squares  for  experimental 
error  and  sampling  error  are  very  nearly  equal  in  all  of  the  analyses  of  variance. 

We  do  not  have  an  explanation  for  this  apparent  phenomenon. 

One  may  note  that  the  magnitudes  of  mean  squares,  particularly  for  sampling 

error,  vary  substantially  from  one  target  area  to  another  in  the  same  table.  This 

results  because  the  A' s  in  the  transformations  as  exhibited  in  Tables  III (a)  and 

2 

III (b)  vary  over  the  target  areas.  The  mean  squares  vary  inversely  as  A  as  an 
approximation  and  this  largely  accounts  for  the  variation  in  sampling  error  mean 
squares  observed. 

We  turn  to  discussion  of  the  analysis  of  covariance  in  Tables  A-4(a),  A-4(a) 
Continued,  A-4(b),  Part  1,  and  A-4(b),  Part  2.  No  significances  for  seeding  result 
for  the  ground- seeded  data  in  Table  A-4(a)  for  any  of  the  five  target  areas  or  any 
of  the  four  models.  The  covariates  do  substantially  reduce  sampling  error  mean 
squares  and  experimental  error  mean  squares  and  increase  the  seeding  or  treatment 
mean  squares,  F-ratios  becoming  considerably  larger  than  in  corresponding  analyses 
of  variance.  Control  area  precipitation  is  the  most  effective  covariate,  followed 
by  duration  or  band  passage  time,  while  the  basic  six  cloud-physics  covariates  are 
only  moderately  effective.  The  effectiveness  of  covariate  sets  is  assessed  by  exam¬ 
ination  of  F-ratios  formed  by  the  ratios  of  covariates  mean  squares  to  corresponding 
sampling  error  mean  squares.  The  F-ratios  for  seeding  are  the  ratios  of  seeding 
mean  squares  to  experimental  error  mean  squares. 

Comments  on  the  analyses  of  covariance  in  Table  A-4(a)  Continued  for  ground- 
seeded  data  with  the  extreme  storm  omitted  are  similar  to  those  for  Table  A-4[a). 

The  covariates  have  been  a  little  less  effective,  perhaps  because  they  provided  a 
major  adjustment  when  the  extreme,  not- seeded  storm  was  included.  Borderline, 
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one-sided  significances  for  seeding  are  achieved  for  target  areas  (i),  (iii),  and 
(iv)  for  model  (3)  with  the  six  basic  covariates  and  duration.  The  possible 
effect  may  well  be  in  target  area  (iii),  the  downwind  area  included  in  both  (i) 
and  (iv).  The  apparent  significances  should  be  discounted  in  analogy  with  the 
problems  with  k^  and  in  the  corresponding  analyses  of  variance. 

Only  covariate  models  (2)  and  (3)  could  be  considered  for  the  aerial-seeded 

data,  since  X  was  not  available  because  no  control  area  could  be  used.  From 
c 

Tables  A-4(b),  Part  1,  and  A-4(b),  Part  2,  it  is  seen  that  use  of  the  covariates 
has  reduced  sampling  error  and  experimental  error  mean  squares  but  treatment  or 
seeding  mean  squares  are  also  reduced  so  that  the  general  effect  has  been  to  reduce 
F-ratios  for  seeding  somewhat.  The  covariate,  duration,  plays  no  effective  role 
for  the  aerial-seeded  data.  The  covariate  analyses  here  are  ineffective  with 
results  rather  similar  to  those  of  the  Phase  I  analyses  for  ground-seeded  data. 

No  clear  explanation  of  the  differences  in  the  effectiveness  of  use  of  covar¬ 
iates  for  ground-seeded  and  aerial-seeded  data  in  the  Phase  II  experimentation  is 
apparent.  The  cloud-physics  covariates  were  measured  deeper  in  the  target  area 
relative  to  seeding  location  for  the  aerial -seeded  data.  Examination  of  the  cor¬ 
relations  in  Tables  A-3(a)  and  A-3(b)  shows  that  the  basic  cloud-physics  covariates 
have  somewhat  higher  correlations  with  transformed  precipitation  means  for  the 
aerial-seeded  data  while  the  similar  correlation  for  duration  is  smaller  for  the 
aerial-seeded  data.  For  the  aerial-seeded  data  correlations  between  the  covariates 
and  the  seeding  variable  Z  are  somewhat  higher.  These  results  are  consistent  with 
the  results  of  the  covariance  analyses  but  do  not  provide  the  desired  explanation. 


Tables  A-5(a),  A-S(a)  Continued,  A-5(b),  Part  1,  and  A-5(b),  Part  2,  show  re¬ 
gression  coefficients  for  the  eight  covariates  of  Model  (5)  for  ground-seeded  data 
and  the  seven  covariates  of  Model  (3)  for  aerial-seeded  data,  together  with  the  re¬ 
gression  coefficient  for  seeding,  the  latter  always  being  positive.  Some  compari¬ 
sons  of  these  regression  coefficients  may  be  made  with  similar  tables  of  Bradley, 
Srivastava  and  Lanzdorf  (1979a).  These  tables  essentially  confirm  comments  above 
in  regard  to  the  covariance  analyses.  X^,  500  mb  temperature,  seems  to  play  a  bigger 
role  in  ground-seeded  experimentation  than  in  aerial-seeded  experimentation  -  this 
temperature  was  regarded  as  important  by  the  experimenters.  X^,  stability  class, 
is  the  second  most  important  of  the  six  basic  covariates,  and  seems  to  have  equally 
important  roles  in  both  ground-seeded  and  aerial -seeded  experimentation.  Control  area 
precipitation  is  important,  as  observed  above,  in  ground-seeded  analyses.  The  role 
of  Xj2»  duration,  is  reduced  when  Xc  is  included  in  the  model  and  X12  is  important 
in  the  aerial-seeded  analyses. 
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6.  CONCLUDING  REMARKS 


Phase  II  of  the  Santa  Barbara  Convective  Band  Seeding  Test  Program  was  a 
relatively  small  experiment  in  weather  modification  when  the  inherent  large  variabil¬ 
ity  of  band  precipitation  and  the  two-part  design  are  considered.  The  ground-seeded 
part  of  the  experiment  was  reduced  to  twelve  experimental  units,  one  of  them  at 
least  suspect.  The  aerial -seeded  part  of  the  experiment  consisted  of  30  experimental 
units  when  four  mixed  storms  are  counted  double  and  22  experimental  units  if  they 
are  omitted.  No  seeding  effects  were  significant  in  the  analyses  of  the  ground-seeded 
data  except  in  one  covariance  analysis  for  target  areas  (i) ,  (iii)  and  (iv)  when  a 
high  precipitation,  not-seeded  storm  was  omitted.  No  seeding  effects  were  signifi¬ 
cant  in  the  analyses  of  the  aerial- seeded  data  until  the  four  mixed  storms  were 
omitted  (See  Table  IV(b))  and  covariance  analyses  removed  the  apparent  effect  of 
seeding  even  in  that  situation,  apparently  because  the  covariates  were  correlated 
with  treatment  as  in  the  analyses  of  the  Phase  I  data.  We  can  only  conclude  that 
both  phases  of  the  Santa  Barbara  study  are  suggestive  of  an  effect  of  cloud  seeding, 
but  that  the  effect  has  not  been  sufficiently  well  established  for  general  scientific 
acceptance.  It  is  the  personal  opinion  of  the  principal  investigator  that  the  exper¬ 
imental  results  have  sufficient  promise  to  justify  further  similar  research  if  im¬ 
proved  experimental  design  as  discussed  below  could  be  effected. 

The  Phase  II  experiment  seems  to  have  yielded  data  that  is  particularly  sensi¬ 
tive  to  the  inclusion  or  exclusion  of  certain  storms.  We  have  seen  this  as  the  one 
extreme  storm  is  omitted  from  analyses  of  the  ground-seeded  data  and  the  four  mixed 
storms  are  omitted  from  analyses  of  the  aerial-seeded  data.  We  have  seen  that  the 
use  of  covariates  is  apparently  helpful  with  the  ground-seeded  analyses,  contrary 


to  indications  from  the  Phase  I,  ground-seeded  analyses,  but  not  helpful  with  the 
aerial-seeded  analyses.  We  can  provide  no  explanation  of  this  situation  except  to 
suspect  that  extreme  values  of  the  covariates,  measured  subject  to  substantial  errors 
also,  may  distort  their  effects.  We  see  no  particular  justification  for  omitting 
the  extreme  storm  from  the  ground- seeded  data  and  only  partial  justification  ""or 
omitting  the  four  mixed  storms  from  the  aerial -seeded  data. 

Brillinger,  Jones  and  Tukey  (1978,  pp.  G-8,  G-9)  discuss  exceptional  events  and 
comment  on  both  orographic  rainfall  and  hail  storms.  We  quote  from  the  report: 

•The  only  way  we  can  see  to  tackle  instances  of  such  problems  of  extreme 
difficulty  is  to  eliminate  the  exceptional  character  of  such  events.  The 
only  likely  way  in  which  this  might  be  done  would  be  to  learn  enough  about 
such  storms  as  phenomena  to  be  able,  from  concomitant  measurements,  to  be 
able  to  predict  their  untreated  behaviour  at  least  moderately  well.” 

It  has  been  our  view  throughout  that  the  major  hope  in  the  improved  design  and  analy¬ 
sis  of  future  weather  modification  lies  in  the  identification  and  measurement  of 
improved  concomitant  variables  or  covariates.  The  cloud-physics  covariates  used  in 
the  Santa  Barbara  experiments  were  reasonably  effective  in  reducing  experimental 
errors,  but  they  were  suspect  since  they  were  measured  in  the  area  of  expected 
response  to  cloud  seeding.  Two  other  covariates,  simple  in  conception,  control 
area  precipitation  and  duration  or  band  passage  time,  may  be  even  more  effective 
if  they  can  be  measured  at  pertinent  locations  unaffected  by  cloud  seeding.  It  may 
be  that  extensive  meteorological  research  is  needed  before  really  good  covariates 
are  identified. 
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In  the  Phase  I  experimentation,  with  convective  bands  as  experimental  units, 
criteria  were  established  for  the  selection  of  "seedable"  convective  bands,  accept¬ 
able  experimental  units,  although  difficulties  in  application  of  the  criteria  were 
encountered.  In  the  Phase  II  experimentation,  with  storms  as  experimental  units, 
the  criteria  were  apparently  applied  to  the  first  convective  band  of  the  storm,  but 
it  is  not  clear  if  they  were  applied  to  the  subsequent  bands  in  the  same  storm. 

The  selection  of  appropriate  experimental  units  appears  to  be  of  crucial  importance 
and  it  may  be  that  failure  of  proper  selectivity  led  to  problems  encountered  in 
analyses  of  the  Phase  II  data.  Proper  selection  may  depend  on  proper  concomitant 
variables.  If  it  is  not  possible  to  select  experimental  units  from  appropriate 
criteria  prior  to  inclusion  in  an  experiment,  proper  criteria,  applied  post  facto, 
could  be  developed  to  determine  whether  or  not  the  unit  is  included  or  omitted  from 
the  experiment,  the  decision  to  be  made  by  a  meteorologist  unaware  of  whether  or 
not  the  unit  was  seeded.  Recall  that  the  meteorologist  determining  raingage  pre¬ 
cipitations  was  deliberately  uninformed  as  to  whether  or  not  the  convective  band 
had  been  seeded. 

We  are  inclined  to  think  that  the  changes  to  the  storm  as  the  experimental 
unit  instead  of  the  convective  band  in  the  Phase  II  experiment  was  a  mistake.  The 
change  complicates  analyses  and  reduces  the  available  number  of  experimental  units 
per  season  substantially.  Only  serious  concern  for  a  persistence  effect  of  seeding 
from  a  seeded  band  to  a  subsequent  not-seeded  band  seems  to  justify  the  design 
change. 

We  continue  to  have  some  concerns  about  the  determination  of  precipitation 
attributable  to  a  particular  convective  band  at  a  raingage  and  suspect  that  this 
may  be  an  additional  source  of  experimental  variation.  We  have  not  been  involved 
in  such  determinations  but  the  procedure  has  been  shown  to  us.  This  seems  not  to 
have  given  the  experimenters  concern,  but  further  study  is  indicated. 
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Given  a  substantial  body  of  data,  further  statistical  analyses  are  always  pos¬ 
sible.  Thought  has  been  given  to  analyses  other  than  those  that  we  have  reported. 

But  it  is  our  judgment  that  analyses  completed  have  exhibited  the  essential  results 
or  lack  of  them.  Concerns  about  multiplicity  of  analyses,  often  expressed  relative 
to  weather  modification  experiments,  are  real.  Some  combination  of  analyses  for 
the  Phase  I  and  Phase  II  data  sets  might  have  been  attempted.  But  the  design  changes 
made  for  the  Phase  II  experimentation  made  this  very  difficult  and  any  attempt  seemed 
likely  to  be  unproductive. 

In  summary,  the  Santa  Barbara  studies  seem  to  have  been  done  carefully  within 
the  state  of  the  art  and  with  full  scientific  integrity.  Experimentation  with  winter 
storms  on  the  West  Coast  may  well  hold  more  hope  for  the  clear  demonstration  of  the 

effects  of  cloud  seeding  than  other  types  of  effort  in  weather  modification.  The 
success  of  future  experiments  seems  likely  to  depend  on  the  use  of  better  covariates, 
avoidance  of  possible  seeding  effects  on  the  covariates,  choice  of  good  control  areas, 
better  selectivity  of  acceptable  experimental  units,  the  use  of  convective  bands  as 
experimental  units,  if  possible  persistence  effects  can  be  discounted,  and  the  use  of 
experimental  designs  that  allow  for  storm  effects  as  a  source  of  variation. 
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APPENDIX  TABLES 
TABLE  NOTES 

Appendix  Tables  are  in  two  parts,  Parts  (a)  for  the  ground-seeded  data  and  Parts 
(b)  for  the  aerial -seeded  data  of  Phase  II  of  the  Santa  Barbara  experiment. 

A-l.  Notes  on  Tables  A-l.  Tables  A-l(a)  and  A-l(b)  provide  summaries  of  raingage  pre¬ 
cipitation  data  for  eight  ground  areas  described  in  the  body  of  this  report.  For  each 
convective  band,  the  average  precipitation  in  inches  from  the  operative  raingages  in 
the  designated  area  appears  in  columns  headed  MEAN  and  the  corresponding  variance  among 
observations  at  these  raingages  is  provided  in  columns  headed  VAR.  The  number  of  oper¬ 
ative  raingages  is  given  in  columns  headed  NUM.  Convective  bands  are  listed  by  storm 
number  (columns  labeled  ST)  and  band  number  within  storm  (columns  headed  BD)  given  in 
time  sequence.  Whether  or  not  a  band  was  seeded  is  indicated  in  the  third  column  headed 
TR. 

In  Table  A-l(b),  it  will  be  noted  that  some  storms  have  been  divided  because  the 
treatment  decision  was  not  applied  to  all  bands  within  the  storm  as  intended  in  the  ex¬ 
perimental  plan.  Storms  91  and  92  at  the  top  of  Table  A- 1(b)  are  examples.  In  our 

analyses,  such  storms  were  treated  as  two  storms.  The  first  two  bands  of  storm  91  were 

treated  as  bands  from  a  two-band,  seeded  storm  whereas  the  third  band  of  storm  91  was 
treated  as  a  band  from  a  one -band,  not  seeded  storm. 

In  Table  A-l(a),  storm  4,  band  1  has  extremely  large  means  and  the  means  for  stonn 


8,  band  2  are  also  quite  large.  These  bands  are  discussed  in  the  body  of  the  report. 


(ii) 

A-2.  Notes  on  Tables  A-2.  Gleeson  (1977)  identified  and  tabulated  concomitant  variables 
for  Phase  I  of  the  Santa  Barbara  experiment.  Similar  tabulations  are  provided  in 
Tables  A-2(a)  and  A-2(b)  for  the  ground-seeded  and  aerial-seeded  convective  bands  of 
the  Phase  II  data. 

There  are  only  minor  changes  in  these  new  tables.  Convective  bands  have  both  a 
Band  Number  and  a  Storm- Band  Number,  the  latter  designation  matching  that  in  Tables  A-l. 
Dates  of  bands  and  beginning  and  ending  times  of  band  passage  at  the  seeding  site  have 
been  omitted  as  has  the  station  of  radiosonde  observation.  Other  concomitant  variables 
have  been  recorded  as  before.  Parentheses  in  the  tables  indicate  less  reliable  entries 
because  some  data  were  missing. 

The  concomitant  variables  are  listed  with  brief  explanations.  See  Gleeson  (1977) 
or  Bvadley,  Srivastava  and  Lanzdorf  (1979)  for  additional  details.  The  concomitant 
variables  are: 

Duration:  Time  in  hours  of  band  passage  at  the  seeding  site. 

Treatment :  Seeded  band-S;  not  seeded  band  -  NS. 

Radiosonde :  Number  of  the  radiosonde  used  to  represent  atmospheric  conditions 
associated  with  the  band. 

Time  (PST) :  Time  of  radiosonde  release. 

Mixing  Ratio:  A  measure  of  water-vapor  content  in  the  air  in  grams  of  water  vapor 
per  1000  grams  of  dry  air. 

700-MB  Speed:  Wind  speed  in  knots  of  700-mb  (10,000  ft.)  horizontal  wind. 

Direction,  700-MB  Wind:  Wind  directi o.i,  degrees  east  of  north. 

Mean  Speed,  Wind:  Mean  speed  in  knots  of  horizontal  winds  from  1000  ft.  to  14,000  ft. 


elevation. 


(iii) 


Direction,  Mean  Wind:  Average  wind  direction  of  average  vector  wind,  degrees 
east  of  north. 

500-MB  Teigperature :  Cloud-top  temperature,  decree  Celsius, approximately  19,000  ft. 
Stability  Class:  UL-unstable  with  low  convective  stability  base,  UH-unstable  with 
high  convective  stability  base,  ST-stable,  scored  respectively, 

1,2,3. 

Showalter  Index:  A  measure  of  atmospheric  stability,  large  positive  (negative) 
values  indicate  stability  (instability). 

Stability  Wind:  Speed  in  knots  of  a  theoretical  wind. 

Direction,  Stability  Wind:  Direction  of  the  stability  wind,  degrees  east  of  north. 
Instability  Transport:  A  time  rate  of  horizontal  movement  of  less  stable  air 
measured  in  knots  squared. 

A-3.  Notes  on  Tables  A-3.  Tables  A-3(a)  and  A-3(b)  show  simple  correlation  coeffi¬ 
cients  among  the  concomitant  variables  listed  in  (3.1)  and  Xc,  band  precipitation 
mean  for  response  area  (v),  the  control  area  for  ground-seeded  data,  the  mean  based 
on  the  Bureau  of  Reclamation  Haingage  network.  Detail  on  the  concomitant  variables 
is  given  above  in  Section  A-3.  The  mean  of  transformed  precipitation  responses  for 
Area  (i)  for  a  band  is  designated  z  in  these  tables  and  Z  is  the  indicator  variable 
for  seeding,  Z  »  1  if  a  band  is  seeded  and  Z  =  0  otherwise.  Pairwise  correlation 
coefficients  are  given  for  the  variates,  Xj,...,  X^*  Xc»  z  and  Z.  Some  comments 
on  these  tables  are  given  in  Section  3  of  the  report. 


A-4.  Notes  on  Tables  A-4.  Tables  A-4(a),  A-4(a)  Continued,  A-4(b),  Part  1,  and 
A-4(b),  Part  2,  provide  analyses  of  covariance  for  the  appropriate  models  of  (3-2). 
The  method  of  covariance  analysis  is  discussed  in  Section  4.  Only  mean  squares 
and  pertinent  F-ratios  are  shown  together  with  corresponding  degrees  of  freedom. 
Complete  analysis  of  covariance  tables  may  be  reconstructed  from  the  information 
given.  In  particular  sums  of  squares  may  be  obtained  and  it  can  be  checked  that 
the  total  sums  of  squares  are  equal  for  analyses  within  tables  and  check  with  those 
of  the  corresponding  analysis  of  variance  tables  in  Section  5.  It  can  be  checked 
that  the  effect  of  adding  covariates  is  always  to  reduce  the  sampling  error  sum  of 
squares,  although  not  always  to  reduce  the  sampling  error  mean  square  because  of 
loss  of  degrees  of  freedom.  Seeding  and  experimental  error  stuns  of  squares  always 
add  to  the  Subtotal  sum  of  squares. 

A-5.  Notes  on  Tables  A-5.  Tables  A-5 (a),  A-5(a)  Continued,  A-5(b),  Part  1,  and 
A- 5(b) ,  Part  2,  give  the  regression  coefficients  for  the  covariates  of  Model  (5) 
for  ground- seeded  analyses  and  for  the  covariates  of  Model  (3)  for  aerial-seeded 
analyses,  together  with  the  regression  coefficient  for  seeding  or  treatment.  Inter- 
comparisons  may  be  made  among  these  tables  and  Tables  A-5 (a)  and  A-5 (a)  Continued 
may  be  compared  with  the  corresponding  Table  A-3  for  ground-seeded  data  given  by 
Bradley,  Srivastava  and  Lanzdorf  (1979a). 


TABLE  A- 1 (a) :  PHASE  II  PRECIPITATION  DATA,  GROUND-SEEDED  BANDS 
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TABLE  A-l (a) (Continued)  :  PHASE  II  PRECIPITATION  DATA,  GROUND-SEEDED  BANDS 
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TABLE  A-2(b) (Continued) :  PHASE  II  CONCOMITANT 
OBSERVATIONS,  AERIAL-SEEDED  BANDS 

BAND  NUMBER  AND  STORM-BAND  NUMBER 


VARIABLE 

178 

32-1 

179 

32-2 

180 

32-3 

181 

32-4 

182 

33-1 

DURATION 

1.22 

1.13 

1.27 

1.00 

1.42 

TREATMENT 

NS 

NS 

NS 

NS 

S 

RADIOSONDE 

330 

332 

333 

334 

336 

TIME (PST ) 

1110 

1430 

(1800) 

1600 

2330 

MIXING  RATIO 

3.7 

4.1 

5.6 

5.9 

5.3 

700-MB  SPEED 

42 

36 

47 

*i4 

26 

DIRECTION,  700-MB  WIND 

241 

236 

239 

230 

225 

MEAN  SPEED,  WIND 

23 

21 

28 

33 

17 

DIRECTION.  MEAN  WIND 

218 

223 

217 

213 

217 

500-MB  TEMP. 

-16.9 

-17.9 

-16.3 

-15.9 

-20.7 

STABILITY  CLASS 

UL 

UL 

UL 

ST 

UL 

SHOWALTER  INDEX 

9 

2 

4 

3 

3 

STABILITY  WIND 

41 

40 

56 

66 

33 

DIRECTION,  STAB.  WIND 

217 

213 

208 

208 

212 

INSTABILITY  TRANSPORT 

-267 

-50 

-40 

292 

-14 

TABLE  A-3(a):  Correlation  Coefficients  among  Means  of  Transformed  Precipitations 
Covariates,  and  Seeding,  Target  Area  (i),  Ground-Seeded  Bands 


is  the  Target  Area  Mean  for  the  transformed  data;  covariates  X,Xj,...,Xp  are  defined  in  (3.1); 
is  the  seeding  indicator  variable  with  Z  =  1  if  seeding  and  Z  =  0  if  no  seeding. 


TABLE  A-3(b):  Correlation  Coefficients  among  Means  of  Transformed  Precipitations 
Covariates,  and  Seeding,  Target  Area  (i),  Aerial-Seeded  Bands 


is  the  Target  Area  Mean  for  the  transformed  data;  covariates  X  . .  are  defined  in  (3.1); 

is  the  seeding  indicator  variable  witli  Z  =  1  if  seeding  and  Z  =  0  if  no  seeding. 


TABU;  A-4(a):  Analysis  of  Covariance  Tables  for 
Transformed  Data,  for  the  Various  Target  Areas, 


•Significant  at  0.05  level. 
•Significant  at  0.01  level. 
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One  extreme  precipitation  storm,  Storm  4  with  one  band,  not  seeded,  has  been  omitted. 
‘Significant  at  0.05  level;  for  seeding,  the  test  is  one-sided. 

‘Significant  at  0.01  level. 


TAIW.Ii  A-4(b):  Analysis  of  Covariance  Tables  for  Models  (2)  and  (3), 
for  the  Various  Target  Areas,  Aerial -Seeded  bands.  Part  1  - 


TABLE  A-4(b):  Analysis  of  Covariance  Tables  for  Models  (2)  and  (3),  Transformed  Data 
for  the  Various  Target  Areas,  Aerial-Seeded  Bands,  Part  2  -  Mixed  Storms  Omitted 
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TABLE  A-5(a):  Regression  Coefficients  for  Model  (5)  for  the  Various  Target  Areas, 

Transformed  Data,  Ground-Seeded  Bands 
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TABLE  A-5(a)  Continued:  Regression  Coefficients  for  Model  (5)  for  the  Various 
Target  Areas,  Transformed  Data,  Ground-Seeded  Bands,  Extreme  Storm  Omitted* 
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TABLE  A-5(b):  Regression  Coefficients  for  Model  (3)  for  the  Various  Target  Areas, 
Transformed  Data,  Aerial-Seeded  Bands,  Part  1  -  All  Bands 


TABLE  A-5(b):  Regression  Coefficients  for  Model  (3)  for  the  Various  Target  Areas, 

Transformed  Data,  Aerial -Seeded  Bands,  Part  2  -  Mixed  Storms  Omitted. 
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